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Abstract 
• ^ , 

^ I We study a new mechanism for hadronic helicity flip in high energy hard exclusive re- 

' actions. The mechanism proceeds in the limit of perfect chiral symmetry, namely without 

any need to flip a quark helicity. The fundamental feature of the new mechanism is the 
breaking of rotational symmetry of the hard collision by a scattering plane in processes 
involving independent quark scattering. We show that in the impulse approximation 
there is no evidence for of the helicity violating process as the energy or momentum 
transfer is increased over the region IGeV"^ < < 100 Gey 2. in the asymptotic 
region > 1000 GeV^, a saddle point approximation with doubly logarithmic accuracy 
yields suppression by a fraction of power of Q^. "Chirally-odd" exclusive wave functions 
which carry non-zero orbital angular momentum and yet are leading order in the high 
energy limit, play an important role. 
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1 Introduction 



The theory of hard elastic scattering in Quantum Chromodynamics (QCD) has evolved con- 
siderably over many years of work. Currently there exist two self-consistent perturbative 
descriptions, each with a specific factorization method for separating the hard scattering from 
non-perturbative wave functions. A well-known procedure using the "quark-counting" dia- 
grams has been given by Lepage and Brodsky A consequence, and direct test, of the 
factorization defining this mechanism is the hadron helicity conservation rule 



where the Aj's are the helicities of the participating hadrons in the reaction A + B ^ C + D. 
The fact that this rule is badly violated in almost every case tested suggests two alternatives. 
One possibility, advocated by Isgur and Llewellyn-Smith 0, is that the energy and momentum 
transfer (Q^) in the data is not large enough for the formalism to apply. However, the data 
also shows behavior close to the model's predicted power dependence on Q^, indicating that 
hard scattering of a few pointlike quarks is being observed. The apparent contradiction has 
led to much discussion, and has even caused some authors to suggest that perturbative QCD 
itself might be wrong. 

The second alternative is that another power behaved process causing helicity flip is present. 
In fact the "independent scattering" subprocess, introduced by Landshoff 0], is actually the 
leading process at very high energies [Q. But it has been assumed that hadron helicity con- 
servation would be the same in the independent scattering process as in the quark- counting 
one, since both involve exchange of hard gluons with large Q^. In general, terms proportional 
to a quark mass, rriq, for example, have been understood to cause helicity flip in either model, 
but with amplitude suppressed by a power of m^/ Q relative to the leading term. Such terms 
seem to be quite small and are probably not a believable explanation of the persistent pattern 
of large violations of the helicity conservation rule. 

Here we show that the independent scattering mechanism predicts high-energy helicity 
non- conservation at a much higher rate. The calculations in momentum space are sufficiently 
complicated that this phenomenon has been overlooked for almost twenty years. Adopting a 
transverse position space formalism introduced by Botts and Sterman 0, we show that the 
details rest on non-perturbative wave functions that should be measured rather than calculated. 
These wave functions measure non-zero orbital angular momentum not taken into account by 
short distance expansions. We argue that the novel factorization properties of independent 
scattering processes cannot practically be reduced to the same ingredients used in the quark 
counting scattering. In any case, it is not necessary to flip a quark helicity: the new mechanism 
proceeds unimpeded in the limit of arbitrarily small quark mass and perfect chiral symmetry 
in the hard scattering. 

This paper is organized as follows. In section 2 we review the derivation of helicity conser- 
vation in genuine short distance processes.. In section 3, we present the independent scattering 
mechanism with special emphasis on non-zero orbital angular momentum wave functions. We 
compute the contribution of these components to a helicity conserving reaction in section 4, 
then to a helicity violating reaction in section 5. These two sections contain our most impor- 
tant results at asymptotic and at accessible energies. Section 6 is an experimental outlook. 
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2 Helicity conservation in short distance dominated pro- 



cesses 

First we review the conventional derivation of hadron helicity conservation p[. The quark- 
counting factorization introduces the distribution amplitude ip{Q'^,x) 0, an integral over 
the transverse momentum variables of the wave function for quarks to be found carrying 
momentum fraction x in the hadron. [| For simplicity of presentation we specialize to a single 
pair of quarks, the meson case. Let ipi^Ti^) be a light cone wave function to find quarks 
with relative transverse momentum and light cone momentum fraction x. In terms of the 
Fourier conjugate transverse space variable b^ separating the quarks, then 

Q 

^{Q\x) = Jd^kT^PikT,x) 



Q oo 

= / d^kT f rf'bTe^^-*^- e'^'^^md hrlx) (2) 



In the second line we have expanded the wave function iplhx, x) to exhibit the 5*0 (2) orbital 
angular momentum eigenvalues m: a complete set for the (nearly on-shell) quarks consists 
of the z-axis orbital angular momentum, the energy, and the z-momentum. Now suppose 
the distribution amplitude '^{Q^, x) is assumed to be a good description of a process. Then, 
whatever the angular momentum of the wave function, evaluating the integrals reveals the 
m = element is the sole surviving term in Eq. (^: 

oo 

^{Q^x) =2nJd\hT\ QJi{\ hr I Q)M\ U) (3) 



This shows that use of (p{Q'^,x) imposes two things: as ^ oo the scattering region is both 
"small", since the region < l/Q"^ dominates in the Bessel function, and the scattering 
is "round", i.e. cylindrically symmetric (Fig. 1). In the absence of orbital angular momen- 
tum, the hadron helicity becomes the sum of the quark helicities. The quark helicities being 
conserved at leading order, the total hadron helicities are conserved. The hadron helicity con- 
servation rule (1) therefore represents an exact symmetry of the quark- counting factorization. 
A crucial question is: does this symmetry of the model represent a property of the entire 
perturbative theory. Or, can we simply assume "s-wave" 5*0(2) wave functions to be the 
main contribution as in a non-relativistic picture? 

The answer to both questions is no. In general, quark wave functions themselves are not 
particularly restricted in orbital angular momentum content, even in the high energy limit. For 
example, in the pion (pseudoscalar meson) case the light-cone wave function may be expanded 
on four Dirac tensors allowed by parity symmetry as 

tPap{x, hr, p^) = {AJ'y^ + 5^75 + CJt15 + D^\^, hh^j^f^ > (4) 

^Our convention does not include a logarithmically varying factor of introduced for renormalization 
group analysis in Ref. |Q, Either convention can be used without affecting our argument. 
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where A^^-Dt^ are functions of the hght cone fraction x, the transverse separation and the 
total four momentum of the meson p^. For the moment we do not discuss dependence on 
gauge fixing and a path-ordered exponential of the gauge potential used to make the wave 
function gauge invariant. The D^r-term carries one unit of orbital angular momentum and 
yet scales with the same power of the "big" momentum p+ as the A^r-term, which is s-wave. 
Since the D^^ term has a bx factor, which can be written in terms of bx^i ± «&t,2 this term 
represents one unit of orbital angular mometum. In terms of power counting, then, the m = 
and m 7^ amplitudes can be equally large. We also note that wave functions are not objects 
to be derived in perturbation theory, but instead represent the long-time, non-perturbative 
time evolution proceeding inside a hadron. The non-perturbative Hamiltonian of QCD does 
not conserve spin and orbital angular momentum separately, but instead generates mixing 
between orbital and spin angular momentum. Finally, there is no simple relation between 
"s-wave" non-relativistic models of constituent quarks, and the pointlike quarks resolved in 
large-Q^ collisions, so no statement can reliably be made about quark angular momentum 
content of hadrons. Thus if a non-zero orbital angular momentum component somehow enters 
the hard scattering — and this is a crucial point — then the long-time evolution before or 
after the scattering can convert this angular momentum into the observed hadron spin. It is 
not necessary to flip a quark spin in the hard interaction, because the asymptotic hadron spin 
fails to equal the sum of the quark spins. Such a mechanism is totally consistent with the 
impulse approximation of perturbative QCD. 

The challenge in high energy hadron scattering is therefore to find those large pro- 
cess(es) in which non-zero orbital angular momentum enters, or in other words, to find those 
which are not "round". It turns out that in any treatment relevant to current energies the 
independent scattering process is not "round" but instead "fiat" (Fig. 1). The subprocess 
is highly asymmetric, showing an extreme dependence on the scattering plane. Rather than 
disappearing in the high energy limit, the dynamics of this process increases its asymmetry 
over a range of Q^. 

The origin of the asymmetry is kinematic. Let rj be the direction perpendicular to the 
scattering plane and ^ be a vector in the scattering plane. In the independent scattering 
mechanism, let two (or more) uncorrelated scattering planes be separated at the collision 
point by a transverse out-of-plane distance 6^. The out-going beams of quarks must coincide 
in direction well enough to make hadrons in the final state. Conserving 4-momenta, for each 
pair of scatterings there are three delta functions of large momenta (scaling like p"*"), and one 
delta function of the out-of plane transverse momentum, of the order of 1/ < 6^ >^/^, where 
<> indicates a typical expectation value in the state. On-shell quark-quark scattering is 
dimensionless and scale invariant (up to logarithmic corrections) in QCD. The overall ampli- 
tude therefore scales like the product of these factors, namely like < 6^ >^/^ (Q^)"^^^. The 

power-counting of this argument is quite well-known; what we emphasise here is the role 
of the scattering plane, namely the breaking of rotational symmetry with the out-of-plane 
direction rj. It is as if hadrons "flatten" under impact in the in-plane direction ^, forming a 
cigar-shaped hard scattering region |Q . This fact is very hard to see in covariant perturbation 
theory in momentum space, explaining why it has generally been overlooked. 

The kinematic violation of hadron helicity conservation by independent scattering raises 
several new questions. It is clear that the usual association of leading-twist (short-distance) 
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and large either breaks down or hinges on dehcate dynamical details. Our approach 
will exploit the fact that leading approximations to any kinematically distinct amplitude are 
always perturbatively calculable. For example, higher order corrections of non-leading twist 
type in the distibution amplitude formalism cannot violate the hadron helicity conservation 
symmetry and will not affect our approach. The first non-vanishing contributions to helicity 
violating amplitudes involve extra partons in the short distance formalism. A gluon embedded 
in the hard scattering, for example, could transfer spin to an outgoing hadron. We need not 
consider such processes, because, as discussed later, they are subleading by a power of and 
perturbatively small since such gluons are "hard". It remains to be shown, of course, that 
helicity violation from independent scattering is not suppressed by the same order. That is 
the main technical task of this paper. 



3 Independent scattering: formalism 
3.1 Kinematical analysis 

Botts and Sterman have considered P] the generic "elastic" reaction Mi + M2 M3 + 
M4, where Mj's are light pseudoscalar mesons, at high energy i/i and large angle center 
of mass scattering angle 9 = arccos(l + 2t/s). Their study has shown that the reaction is 
dominated by the 2 independent scatterings of the valence constituents with a kinematical 
configuration depicted in Fig. 2. One has two scattering planes separated at the collision 
point by a transverse distance b. 

To make it clear, let us consider a light cone basis (f «, f •, ^i, i]) attached to each meson Mj, 
and chosen so that, in the center of mass frame where pi = 3 and = cos 6* 3 + sin 6' 1, 

^1=^2 = 71(0 + 3); v[ = V2 = j=^{6-3)- 

6=6 = 1; r] = 2; 

vs = v'^ = ^{0 + sinO i + cos 6 3); = f 4 = ^ (0 — sin 6* 1 — cos 6 3) ; 

6 = 6 = cos 61— sin 6 3. 

In the cm frame, one has Pi = Qvi (Q = Js/2). Each Bethe-Salpeter amplitude Xj 



is assumed to eliminate qq configuration with relative 0{Q) transverse momentum (along 
and rj). The relative minus momentum (along v^) is of 0{M'^ / Q). Then, momentum transfers 
in H or H' are dominated by large invariants built with XjQ-terms and a first approximation is 
to neglect all but these components of the quark or antiquark momenta, in the hard amplitudes 
H and H' . This is the impulse approximation. 

A second observation follows from kinematics; although momentum conservation at the 
hard scattering H relates the internal momentum dependence of the X's, a variation of one 
such momentum fcj in its S,i or v[ direction induces negligible modifications in the three other 
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X's. Consequently, momentum components of fcj along or v[ only appear as a relevant 
variable in the wave function Xi and integrals over these components can be carried out. The 
components along 77, denoted /j, represent transverse momentum out of the scattering plane 
and do not share the same property. Thus, the vertex delta function may be reexpressed as 

a/2 ^ 

5^{ki + k2 - h - ki) = . Y[ 6{xi - Xi)6{li + h - k ~ h), 

I sm u\LJ 

which indicates that the 4 constituents which enter or leave each hard scattering carry the 
same light cone fraction {xi = x 01 1 — x). Introducing the "out of plane" impact parameter 
b through 

r+00 



one may write the amplitude of the process P as 

\/20 /■! r+°° 
K^^t) = ] ^ / dx{27c)^Hi{xQv})H'i{xQv}) / dbl[Vo.^pX^,b,Qvi). (5) 

ZTT Smt^ JO J-00 1 

where "P is the amplitude 

V^f^ix, 6, Qv) = e'^^y- < 0| T {qM 1^) KQv) > _ , (6) 



27r 



y=y-v'+hri 



with Dirac indices a[3. Color indices are suppressed in Eq. (j^) and sums over repeated indices 
are understood. We consider unflavored quarks. The implementation of flavor is straightfor- 
ward by setting to some of the graphs we are going to consider. H and H' are Feynman 
amplitudes (a sum over allowed diagrams is understood) at the level of quarks computed with 
standard perturbative QCD vertices and internal propagators; they consist, at lowest order in 
the coupling constant, of one-gluon exchange or qq annihilation for each quark pair. Note that 
V and H are not individually gauge invariant. If the short-distance 6—^0 limit is assumed, 
then the four general wave functions in Eq. (^) can be reduced to one A^r-term, by taking 
the trace Vq of Vap with Qf'~i%p. For the case of helicity conserving amplitudes, this selection 
was shown to be self-consistent by Botts and Sterman]^. Then, the zeroth moment of Vq is 
related to the decay constant of the corresponding meson 



/ dxVoix, b = 0,p = Qv) = fM 
Jo 



where, e.g. for the pion, = 133MeV. This zero-distance quantity contains no information 
on the interesting dependence on the transverse variable b. 

Here we are concerned with the leading order description of helicity violating terms. Thus, 
we will consider A^-type and Z^Tr^type amplitudes on an equal footing, and make no a priori 
assumption that the region & — > dominates. 



3.2 Gauge Invariance 

The development so far has been sufficient to isolate the kinematic region of interest, which as 
we have already noted is characterized by finite separation between the participating quarks 
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in the out-of scattering plane direction. The amphtude is thus a strong function of the spa- 
tial dependence of the wave function. The Bethe Salpeter wave function is a bilocal matrix 
element and is not gauge invariant. However, we will now discuss how gauge invariance of the 
description can be obtained. 

The key is in how the perturbation theory is re-arranged. In the Sterman and Botts 
factorization certain "soft" corrections are put into the wave functions, leaving other parts of 
Feynman diagrams to go into the hard scattering kernel. In dressing the wave function in this 
way, it is no longer a quark correlation (the Bethe Salpeter wave function), but the matrix 
elements of operators determined by the types of diagrams put in. The operators chosen in ^ 
are path ordered exponentials (poe's), shown by Collins and Soper to be the generators of 
eikonal approximations to the gluon attachments. The poe's are gauge covariant, leading to 
a gauge invariant amplitude. 

This is partly forced by physics, and partly a convention. As a convention for the per- 
turbation theory, subsequent diagrams must be evaluated with subtractions to avoid double- 
counting. More generally, any operator functional of the A fields which transforms properly 
could serve in place of the poe's, and creating a different subtraction procedure. Let us extract 
what we can that is independent of convention. 

Let the operator in the definition of the wave functions be called U{A] x); we will call it a 
gauge-dressing operator. Under a gauge transformation at the position x, we require U{A; x) 
to transform like an antiquark. Then products such as U{A] x) iplx) are gauge invariant. That 
is, we have gauge invariant matrix elements to find a dressed quark 

{0\T{U{A; y)^{y) U^{A; x)ij{x))\7r{p)) 

It is obvious that this requirement does not determine U {A; x) uniquely, because one could 
always attach a factor which is gauge invariant without changing the gauge transformation 
properties. The particular choice of what to attach is a prescription, i.e. a definition of what 
parts of the amplitude will be put in the wave function and what in the hard scattering, and 
it cannot be determined by gauge invariance alone. However, due to gauge transformations 
one must attach some kind of gauge-dressing operator to have a well defined matrix elements. 

3.3 Path-Independent Dressing 

Although the standard way to do gauge dressing with the poe is path dependent, no path 
dependence generally need be associated with U{A] x) and in particular the observable process 
does not determine or favor any path. This important point can be seen with an elementary 
example from QED, where the U{1) gauge invariance is easier to control. The straightforward 
QED analogue of our process involves equal time (not light cone time) correlation functions 
in the gauge ^4° = 0. This gauge choice eliminates a mode, but there still remains a lack of 
definition of the coordinates due to time independent gauge transformations 9{x) 

A{x) A'{x) = A{x) + Ve{x), d^e = 

This produces a change in the longitudinal modes. These modes are sometimes also called 
unphysical, a very unfortunate choice of terminology. In free space and in the absence of 
couphng a gauge theory has two transverse degrees of freedom and the third would be called 
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unphysical. However we are interested in the case that matter fields exist (and the non-Abehan 
couphng is turned on) in which case the third mode is real, but special inasmuch as being 
constrained in terms of the other variables by gauge invariance. To see this, note that we can 
decompose into transverse and longitudinal parts, 

A = At + Al 
= At + V0 

(x) (transforming part) 



At — A — V0 (invariant part). (7) 

Since it is there for gauge transformations, the longitudinal part (f) is not free to be varied 
in arbitrary dynamical ways, but must accompany the matter field in a prescribed, unique 
functional. At time t — d this operator is 

U{A; x) = e*»'^(^) (8) 

from which one can verify that under a gauge transformation, 

0(x) (i){x) + e{x) 

iPix) e-*»^(^V(^) 
U{A;x) e^^^(^)f/(A;,x) 
U{A;x)ilj{x) U{A;x)ilj{x) (9) 

We will call this the Coulomb dressing because it creates a classical Coulomb field around 
the matter particle, as the reader can check by calculation. Since the Hamiltonian commutes 
with the gauge transformation operator once A^ = has been set, the time evolution will 
maintain the invariance of the combination U{A]x)il){x). However, as noted already, one is 
not forced to accept this as the unique answer, but can opt for [7(^4; a;) /(At), which will 
transform in the same way for any /(At). 

The reader may still be curious to know the relation between Coulomb dressing and the 
poe approach. This can be very simply exhibited by noting that 

^(^)e«fl<^(2/)e-i90W^(a.) = ^jj(y)e'9 1" '^-■^ ^jj(x) (10) 
This expression is still path independent. This is the choice /(At) = 1. A different choice 

is 

/(At) = e'^-/'>"'^^(^) (11) 

in which case we have 

U{A;y)U^iA;x)f{AT) = e'^T.'^^M-) (12) 

which is the standard path dependent poe. Both procedures are equally acceptable, as far 
as satisfying gauge invariance, but the path ordered exponential can create a line of physical 
transverse gauge field particles between the charged matter fields, depending on how it is 
oriented. Such a line of gluons, existing only along the chosen part, can be interpreted as 
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an arbitrary model for the transverse gauge field inside the state of interest. Similarly, if one 
boosts the Coulomb dressed definition, the boost also creates a blast of transverse gauge fields 
as seen by the "equivalent photon" approximation. 

In perturbation theory, the lowest order approximation to non-Abelian dressing is the 
Abelian case. It is possible in the non-Abelian theory to write down expressions analogous 
to the Coulomb dressing but care must be used to keep track of the color indices. It is 
equally valid to use poe's, which definitely transform properly, as building blocks to generate 
an infinite number of different ways to dress the quarks. The different choices are not relevant 
for a leading order calculation to which we restrict this study. 

3.4 Factorization 

The next crucial step is to elaborate a factorized form for the amplitude, whose prototype 
is Eq. (^, regarding radiative corrections. Generalizing the results of |^ to the case of the 
helicity-violating Dirac projections, a leading approximation to the soft region rearranges 
these corrections to obtain the following expression 

= ,1 .^^1 / dxH{{xQv})H'{{xQv}) dbU{x,b,Q)m'\x,b;l/\b\)) , (13) 

Ztt\ smt)\ Jo J-i/A ^ ' 

where B. and H' are evaluated at respective scales xQ and xQ which are assumed to be 
large. Large logarithmic corrections to the process, with the coexistence of the two scales 
Q and 1/6, are resummed in f/, in such a way that is now a soft object: it does not 
include perturbative corrections harder than 1/6. It is the non perturbative object necessary 
to connect short and long range physics, both present in hard hadronic processes. 

Expressions for the above quantities will be given in the phenomenological study, but let us 
specify here some of their qualitative features. When b is small, smaller or of the order of l/Q, 
radiative corrections are also small and must be considered as perturbative corrections to the 
hard amplitudes H or if'; in this regime, V^^^ is closely related to the distribution amplitude 
ip{x] fi) evaluated at the renormalization scale /i = 1/6 0. Depending on the definition, the 
distribution amplitude also includes resummed logarithmic corrections from hadronic scale up 
to /i [|I], 0- When 6 is large, this replacement is not legitimate but since f/(6) is a strong 
suppression factor in this region the exact value of P^"-* is irrelevant. 

Endpoints in the x integral where hard subprocesses would indeed become soft may look 
problematic. Both the distribution amplitude and wave function approaches used here become 
self-consistent because of the end-point zeroes, e.g. ip{x — > 0) ~ x'^, where k > which should 
occur independently of spin projection. 

4 Contribution from non-zero orbital momentum com- 
ponents of the wave function: tttt nn 

Before analyzing helicity violating processes let us examine the leading contribution from the 
various components in Eq. to a standard helicity conserving reaction such as vrvr — * vrvr. 
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4.1 Computation of hard amplitudes 

In their study, Botts and Sterman were interested in identifying the asymptotic behavior of 
the amphtude A{s,t) (s +00, ^ fixed). Asymptotically, the Sudakov mechanism contained 
in U{b) results in a suppression of large b region in the integral of Eq. (|1^). In this limit one 
can forget about tensorial components of Va/3{x, b, Qv) oc {■ ■ - p- ■ ■)a/3 and only the projection 
of Pa/six, b,Qv) onto the particular tensor |75^|q,/3 survives. 

In the intermediate regime, configurations of the qq pair sitting in a light meson with 
transverse separation smaller than the meson charge radius are not strongly affected by the 



Sudakov mechanism |T0|. As anticipated in section 2, m = 1 components of the wave function 
which form large invariants in H (as large as the s-wave term) may give sizeable corrections to 
the interaction amplitude between pure s-wave states. For contributions with leading power 
behavior in the pseudoscalar case, we must keep the tensorial decomposition 



Va(^{x,b,Qv) = ^-f5{Vo{x,b,Qv)i, + Vi{x,b,Qv)[^,P]}^f,. (14) 



We now explore the calculation with this assumption. To begin with, one forms the 
projection, denoted as t and t', of the hard amplitudes H and H' on the relevant various 
Dirac tensors coming from the wavefunction decomposition. We follow Botts and Sterman's 
classification of graphs, with three fermionic flows for H and two gluonic channels each (see 
Fig. 3), 

/ H{M^Mj^ M^M^) gluon c(^^^^ c{^^^^ 

1 12 > 34 M, S ^a\a2^a3ai ^a\aii^a2a3 

2 12 ^ 34 t, S ^a\a2^aza4, ^a\az^a2a4, 

3 12 > 34 t, 11 5a\ai^a2a3 ^0103^0204 

Color flow in this problem is simplified by noting that one gluon exchange between two quark 
lines gives a color tensor 



which we may reexpress with the color tensors listed above, in the form 

2 2 _ 

1=1 1=1 

with ^ ^ 

Ci = €2 = -, C2 = Ci = —-. 

I 

With this notation, one finds for the hard amplitude Oq, containing no b factor, 

su 

.(2) ^ ^ s'^ + . ,^ ^ . ^ ^ .u^ . ^ ^ + 



{titjY," = CiCj—^ + {CjCj + CiCj)- + CjCj—^; (15) 

st 

{titj)f = CiC/^^^ + {CiCj + CiCj)'^ + Cic/^.'^ 

V tu u 



j- 



,2 
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times an overall common factor 



- 



,6/ x'^x'^s'^^ 

where g is the QCD coupling constant which appears in Feynman rules. We have already 
indicated that the whole amplitude Eq. can be properly defined regarding renormalization 
and factorization, so that g'^ stands for {4TT)'^as{xQ)as{xQ). 

There is no or terms, due to the odd number of 7 matrices; this is a consequence of 
chiral symmetry. The second term is therefore a hard amplitude 02, containing 6^, which is 
found to be 

{tjtjYi^ = CjCj^ + iCjCj + CjCjf^^ + CjCj^-^, (16) 



with a common factor 

(here 6 is a distance so that 6^ > 0) and the third one a 6^ hard amplitude, a^, which is the 
combination of Eq. (|1^) , but with an overall factor 

![y 5iv^4_ 

6/ x^x^s^ 

A check of the above expressions or a possibility to reduce the number of graphs one has 
to compute is provided by symmetries under meson exchange; starting from the expression 
one gets with two w-gluon exchange and fermionic fiow / = 1, see Fig. 3, which we label uul, 
one can generate 

Exchange channel kinematic color 
2 <-i> 4 ssl u ^ s C ^ C 

3^4 tt2 u^t C 

2^3 mm3 s t C C . 

The reader will easily find the channels obtained from another starting point, say msI, and 
the combination of exchanges needed to determine graphs which do not appear in the above 
array, thus completing the whole amplitude. 



4.2 Asymptotic behavior 



We are now ready to evaluate the integral defined in Eq. (13) with the above hard amplitudes. 
In a first step, we approximate the Sudakov factor by its dominant expression at large Q H] 



U{x, b, Q) ^ exp 



xQ 

-cln ~^(~ lnu(xQ, 6) — 1 + u{xQ, b)) 



exp [x ^ x] , 



:i7) 
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with 



4 2 

and c = 4 — = ?)2/2'l for nj = 3. We have introduced the variable u{x,b) which 

3 11 — 2,nf /3 

turns out to be the relevant one to describe the Sudakov unsuppressed region in the (6, Q)- 
plane: for u{xQ, b) = 1 there is no suppression from the first exponential in Eq. |l3; as soon 
as u departs from unity, one gets rapidly a strong suppression due to the large InxQ factor. 
Eq. (0) is likely to be valid for any Dirac projection because the leading logs factor away 
independent of spin. In this approximation, U is a. scalar in color space and one easily performs 
the color traces, with 



Cl{a}Cj{b} n ^ci^b, = 3 ^ ^ 3 j , 

1/9 -3 \ ^ ^ _ 1 / -3 9 \ 
^^^-^ - 36 I, -3 1 j ' ^^^^ - 36 i 1 -3 j ' 

^^^^ = ^(-3 9^)' 
and obtains the following hard amplitudes, labeled by the power of b entering: 

Try 256^^ s\s^ - 3tu) + t^u^s^ - tu) _ 
6/ x^x^s^ sH'^u^ ' 

Try 2048^^ - 2,tu) - t^u^ 
6/ x^x^s^ s^t'^v? ' 

TT^ 2048^^ s\s^ - 3tu) + eu\s^-tu) _ 



Oo = 

«2 = - ) 

04 = 6^ 



which at 90° is 



6/ x^x^s^ s^t'^ 



ao = 19 - ^ ; 02 = 1206 - 



67 x^x^s^' v6/ x^x^s^' 

\6/ x^x^s^ 

For the remaining soft wave function Vq, one can adopt two different models. First, one 
can approximate Vq to its value at small b (see and the argument in section 3.1) 

Voix,b;iJ,) ^ ip{x;iJ.), (19) 

with 

Vas{x) = 6/^x(l-x), (20) 

ip,,{x;fi^ 500MeV) = 5(2x - l)Vas(a;), (21) 

which are standard choices for pion distributions. The asymptotic form Eq. (^) is derived 
in [|ll, 0]; it has no evolution with fj, and is indeed the limit as /z — oo of all distributions. 
Because evolution with fi is slow, the effective distribution at non asymptotic regime may 
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be very different from ipas- Chernyak and Zhitnisky |Tl[] have built from QCD sum rules 
the above cz-form Eq. (|2T|) , which evolves with the scale, but with quite small effects on any 
computation when one does not examine a huge interval of energy. We will ignore these effects 
in the following. 

The existence of the additional component Vi complicates the problem. Next to nothing 
is known about it, but a reasonable ansatz is to adopt a form similar to Vq. One notices that 
Vq has a mass dimension which appears to be set by In the case of Pi, the dimension is a 
squared mass and we do not know what the normalization constant has to be. We will adjust 
the normalization in the way described below, and assume the same x-dependence for Vq as 
Pi. 

For each hard part Oj, freezing the coupling at the scale Q to simplify the study, we get 
through Eq. (|T3p a value Ai{Q, 90°) and perform the ratios Ri = Ai/Aq. 
To get the asymptotic behavior, we analytically evaluate the 6-integral 

rf66"t/(6,x,Q) 



with a saddle point approximation; using the change of variable u = — In 6/ In y/xxQ. One 

2c 
2c+n+l 



has a maximum of the integrand at uq = „ . and finds 



^ dbb''U{b,x,Q) ^uo\j^^-^{xxQ^y'''''°. (22) 
Performing the x-integration 



i{ip,n) = f\xip{xY {x{i - xyy^'^"^"^^^^ , 







an asymptotic expression is found for the ratio 



- 2c + n + l [a) 7{^y ^^^^ 



from which one deduces 

i?2 oc g-^-^°, i?4 oc g-^-^^ 

Roughly speaking, each power of b is suppressed by 1/ y/Q (not l/Q). 



4.3 Intermediate behavior 

At accessible energies, we expect deviations from the result given in Eq. (^3|) . This comes 
first from an eventual failure of the saddle point approximation with too small Ing. We 
thus have numerically evaluated the amplitude Eq. (|l^) with U given by Eq. (|l^). Results 
for our computation of the ratio of amplitudes are displayed in Fig. 4 with cz distribution 
(solid line). We get similar results for the asymptotic distribution amplitude (Fig. 5). To fix 
the normalization, we choose here and in the following to set arbitrarily the ratios to 1 at 
y/s = 2GeV. We observe that R2 decreases by a factor around 7 from ^/s = 2GeV to 20GeV. 
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This is a much milder suppression than the naively expected factor (dotted curve). i?4 
drops more drastically by a factor around 20 in the same energy interval. A numerical study 
shows that logarithmic corrections ignored within the saddle point approximation are not 
negligible in the accessible range of energies. At larger values of the energy, y/s > 20GeV, the 
approximated result of Eq. ( ^5] ) becomes accurate. 

Secondly, Eq. (^) should be supplemented with non leading terms which are known in 
the s-wave case Q but presumably different in the p-wave case. The neglected logs in the 
expression of U Eq. ([T7|) may modify soomewhat the ratio over some intermediate range of 
energy. To modelize such an effect, we add a simple x, 6-independent term like 

exp In In ^ 

in the expression of R2 with K some constant. The ratio modified by such a factor is shown 
in Fig. 4 as a shaded area limited by the curves corresponding to = 1 and K = — 1; 
this measures in some way the theoretical uncertainty on p-wave contribution suppression. 
Further theoretical progress in the computation of these In In Q terms will allow to get rid of 
this uncertainty but we do not tackle this task here. 

A third effect may come from the intrinsic transverse dependence of the wave function. 
While the replacement of V^^\x,b; 1/6) by Lp{x; 1/6) discussed in section 3.1 is reasonable at 
large Q (> lOGeV), it is more questionnable at intermediate values where long range physics 
may be accounted for by including some intrinsic 6-dependence |Ty] as 



Vo{x, h) = 4nAf^{x) exp (-^ - ^) . (24) 

with parameters = .096, = .88GeV^ and A/" = 1.68 for the asymptotic distribution 
amplitude. The results are modified as depicted in Fig. 5 where the curves from this wave 
function and from the asymptotic distribution are shown for comparison. Furthermore, we 
notice that the phenomenology may also be modified considerably by allowing variations of 
different distribution form, that is yjas for ^0 aiid ipcz for Vi or the reverse choice. 

These three effects show that the power-like decrement Eq. ( |23D of the ratio is diluted at 
intermediate energies and consequently the amplitudes Aq and A2 are likely to compete over 
a rather large interval of s, say IGeV^-lOOGeV^. 



5 Helicity violating processes: nn pp (hs + /14 7^ 0) 
5.1 Double helicity flip 

To begin with, let us find the possible tensorial decomposition for the quark antiquark wave- 
functions of the p meson. One specifies a one particle state by the momentum pi {i = 3,4), 
which with the notation of section 3, one can write: Pi = Qvi + f^^i (m is the p mass which 
we do not neglect for the moment), and by the helicity hi G {1,0, —1}. The 3 helicity states 
are described in a covariant way with the help of the 3 vectors e'^{pi), satisfying Sh-Pi = and 
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efi-£h' — ~^hh'- For pi^ps), {'P3,C,3,f]) is a direct trihedral, and we choose 

(^ = +1) £r = -^i^3 + ^v), 

ih^ 0) .0 - £.3-^.^, (25) 

(^ = -1) eL = 



V2' 



A consistent choice for p{p4) is 



(^ = +1) £r = -^i^s-iv), 

{h= 0) .0 = (26) 

(^ = -1) eL = -^{^3 + ^v)■ 

Thus, demanding a parity — state, the Bethe-Salpeter amplitude has the most general 
Dirac-matrices expansion 

r{x,b;p,h)= (27) 
Vo Ih, f\ + ^0 4 + ^1 ^h-hi) + V[ Sh-h + Vi [^hj] H'P'i l^h, f\+V2 Sh-h lf\+V'^ Sh-h l/>. 

One can then extract the relevant components for the study of independent scattering 
processes. That is, we isolate dominant high-energy tensors, those which contain one power 
of the large scale Q and we restrict to fej. {i.e. bq^ in the notation of section 3.1). We get for 
a longitudinally polarized p 



V{x, br; p,h = 0) = ^(v'J + V[ M) , (2^ 



Q 
m 

and for a transversally polarized one 

V{x, 6t; p,\h\ = l) = Q (Vo [4, i,] + Vi Sh-hr i' + Vi [4, i^] ^ + V2 eh-hr [?^, h]) ■ (29) 

Before going through the computation of an helicity violating process, let us exhibit what 
we call the double-flip rule. It follows from the property of each dominant component that 
each power of b changes the chirality (which is + if the component anticommutes with 75 and 
— if it commutes). With vector gluon couplings, it follows that if the process is allowed by 
hadron helicity conservation, chiral symmetry requires that the 6-components occur in pairs.0 

^The power suppression for amplitudes with extra transverse gluon in initial or final Fock states follows 
from the same consideration. By power counting one extra gluon embedded into one of the 2 hard scatterings 
is down by 1/Q with respect to the short distance contribution from valence states. However the chirality 
of the component is the same as the p-wave we are considering (that is opposite to the s-wave) . The above 
argument then applies to ensure that when p-waves and one-gluon components are likely to compete, at least 
2 such components are necessary: from section 4, one counts Q~°^^^ for one p-wave and for one extra 
gluon; this shows that considering one extra gluon in the Fock state is indeed subdominant with respect to 
considering a p-wave state. 
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We have already encountered this property in the tttt elastic scattering. Furthermore, for 
an helicity 0-state, in both pseudoscalar and vector case, the chirality of the s-wave term is 
positive, whereas for an helicity odd state, it is negative. Then with chiral symmetry and 
because the total power of b has to be even, one concludes that the number of helicity-0 
states has to be even too. For a 2 to 2 process, only those processes that violate the helicity 
conservation rule by two units are expected to be non-zero in the present framework and hence 
important at large energy. 

An example of interaction which satisfies neither the helicity conservation rule nor the 
double-flip rule is 

Tin PqPt] 

such an amplitude is thus vanishing by power counting. 

We then turn to the study of an helicity violating process which, however, satisfies the 
double-flip rule, namely tttt PrPr. 



5.2 TTTT PrPh 

Let us first verify of the vanishing of the hard amplitude using the s-wave components of the 
external mesons. For this purpose, it is instructive to examine the connection between quark 
helicities and the s-wave Dirac tensors we have used until now. This is accomplished in the 
following way. One considers the free massless spinors of a quark and an antiquark moving in 
the same direction, the quark having a momentum xp and the antiquark a momentum xp, so 
that the compound system has a momentum p. Then one constructs the four possible helicity 
states of the system with solutions of the Dirac equation and finds 

(vr) -^{uc{xp,])vp{xp,i)-Ua{xp,[)vp{xp,])) = 

(Po) -^{Uo,{xp,'l)V(3{xp,i)+Uaixp,i)vp{xp,'\)) = 

(Pr) Uaixp, Dv/^ixp, t) = 

(Pl) Ua{xp, i)v,3{xp, i) = \l'^h'i>Up^ (30) 

The helicity conservation rule is then easily verified with these combinations of spinor when 
one chooses the chiral representation |T^. In this representation, the 2 diagonal blocks of each 
7^ are equal to the null 2x2 matrix. A Feynman-graph fermion line, with vector (or axial) 
couplings and massless propagator, is an even number of 7-matrices between two spinors: 

,lj'\p',h')^0l^^^^--4nl^-+^'lp{p, h), 

Inclusions of 75, which is diagonal in this representation, do not modify this property. Then, 
since to order m/p the chirality of a spinor corresponds to the helicity of the state and, in this 
representation, chiral eigenstate spinors have either their 2 first components or their two last 
equal to 0, states of different helicity always give a null product. 
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This property is algebraic and therefore independent of the representation chosen. How- 
ever, it is more difficult to observe its effects in the trace formalism described in section 2. 
Let us examine how it works in the case of the reaction considered. Among the twelve graphs 
discussed in section 2, each tt and uu graph is because they contain traces over the product 
of an odd number of matrices. For the 8 remaining graphs, the sequence 

occurs (the anticommutation to the left of every 75 is understood and does not modify the 
reasoning). The product of an odd number of 7^^ being a linear combination of 7^^ and 757^^, 
one is left with the evaluation of 

i^ni^l^hf ■■■ and i^r f 757^ ^if 

where = v'^ and e^R = e^i have been used and the index 3 dropped. One can decompose 
each onto (^, ^r, ^l), which are such that their square is and their anticommutation 
rules are 

{^,h} = {^,h} = {f,^R} = {f,h} = o, 

to conclude that all graphs effectively vanish. 

Adding the contribution from one p-wave function does not give any contribution (even 
before the integration over b), because the total number of 7^ is necessarily odd. The first non 
zero term is therefore a 6^ hard amplitude M2 and the computation, within the approximation 
of section 4.2, leads to 



TT 



+ 4^^ol^iP^3p + 4(.^ - 3t« + ^ - 2'^)VlVl\ . (31) 



This amplitude M2 has to be supplemented with U (Eq. {\L7\)) and integrated over b and x. 

Even though this combination involves several unknown objects, one notices that the an- 
gular dependence varies from one component to another and is rather different from the one 
obtained in tttt elastic scattering. Therefore, it may be possible to analyze the contribution to 
helicity violation processes from different wave functions and use this information to deduce 
properties of the wave functions. 

The numerical study of section 4 can then be used to understand the energy dependence 
of double helicity violating processes at accessible energies. As explicitly shown in Fig. 4a and 
5, the naive 1/Q^ factor is replaced by a milder suppression, and this is primarily due to the 
specificity of the independent scattering mechanism supplemented by Sudakov effects. Even 
at very large energies the Q~^'^^ ratio of Eq. (|23|) looks quite weak a suppression. 
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6 Realistic processes and experimental outlook 



Studying meson-meson scattering is an unrealistic simplifying assumption. Including baryons 
is a necessary but quite intricate further step, because of the high number of Feynman dia- 
grams and internal degrees of freedom to be integrated over when more valence quarks are 
involved. We can however still draw some conclusions from our analysis, leaving to future 
work a comprehensive study. The mechanism we have explored occurs in several experimen- 
tally accessible circumstances. Indeed, there is a host of reactions involving hadronic helicity 
violation from which we could learn about the interface of perturbative and non-perturbative 
QCD. 

The helicity density matrix of the p meson produced in vrp — > pp at 90° is a nice measure 
of helicity violating components. Experimental data fl^ yield pi_i = 0.32 ± 0.10, at s = 
20.8GeV^, 6'cM = 90°, for the non-diagonal helicity violating matrix element. Without entering 
a detailed phenomenological analysis, we may use the results of section 5 through the following 
line of reasoning. Assuming that the presence of the third valence quark, which is not subject 
to a third independent scattering, does not alter much the results, one may view pi_i as 
coming from the interference of an helicity conserving amplitude like vrvr PlPr with a 
double helicity flip amplitude like vrvr — prPr. We then get a mild energy dependence of 
this matrix element i.e. Q^^'^^ (Eq. (|23|) ) at asymptotic energies or as shown in Fig. 4a and 
Fig. 5 at accessible energies. This is at variance with the picture emerging from the diquark 



model [|14|. Measuring the energy dependence of this effect is thus highly interesting. 

The most well-known example of hadron helicity violation occurs in pp pp scattering 
Our demonstration of helicity violating contributions to meson-meson scattering has a bearing 
on this, because generalized meson scattering is embedded in the diagrams for proton scat- 
tering. Without needing to make any dynamical assumption of "diquarks" , the perturbative 
QCD diagrams for pp — > pp scattering contain numerous diquark regions, convolved with scat- 
tering of an extra quark. There is no known selection rule which would prevent the scattering 
of such a subprocess with helicity flip from causing helicity flip in pp elastic scattering. This 
does not exhaust the possibilities, because there are other channels of momentum flow and 
color combinatorics which might have different interpretations. The data for pp pp also 
reveals large oscillations about power-law behavior, a second piece of evidence that the short 
distance picture is inadequate. Elsewhere [16| we have identified these oscillations as a sign 
of independent scattering. Given the theoretical and experimental evidence, we therefore 
find no evidence that hadronic helicity conservation is a feature of perturbative QCD, and we 
believe that independent scattering is a main contender in explaining the observations. 

Since reactions of baryons are extremely complicated, and next to nothing is known about 
the various wave functions in the proton, a productive approach to the question is to ask 
for experimental circumstances in which the general mechanism we have outlined could be 
tested without requiring too much detail. We believe that progress here will come from using 
nuclear targets, and studying the phenomenon of color transparency in hard (as opposed to 
diffractive) reactions. This program has been outlined elsewhere [|^; it suffices to mention 
here that suppression of large 6^ regions is expected in reactions of large nuclei. It follows that 
helicity conservation should be obtained in the same circumstances. Thus the mechanism we 
have outlined is testable. It is time to go beyond the short-distance physics used to establish 
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QCD as the correct theory, and we beheve that a multitude of phenomena involving spin, 
color transparency, and detailed hadron structure, will play a major role in the future. 
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Figure 1: Coordinate space pictures of meson-meson independent scattering, (a) Trajec- 
tories of quarks. In the scattering plane quarks approach each other within a distance of order 
l/Qi while the transverse separation of scattering planes (indicated in perspective view) is 
larger and set by the wave function, (b) Contour map of real part of the product of m = 
leading Sudakov wave functions (from |]^) which give a model for the integrand of Eq. (|^). 
No other soft wave function has been introduced, (c) Same wave function multiplied by a 
polynomial representing m = 2. Contours in (b,c) show flattening in the ^-direction of the 
effective wave function imposed by a Gaussian with 6| < l/Q^ at =2 GeV^; higher Q"^ 
increases the flattening. Units of ^ and are fm. 

Figure 2: Kinematics of the independent scattering mechanism. 

Figure 3: Feynman graphs for the lowest order hard amplitude if; for H' reverse the 
arrows. 

Figure 4: The energy dependence of the i?2 (a) and of the (b) ratios with cz distribution 
amplitude 0] (thick lines) and their naive behaviors (thin lines), respectively l/Q^ and l/Q^. 
The shaded area in (a) indicates uncertainty from neglected logs. The ratios are normalized 
to 1 at =2GeV. 

Figure 5: The energy dependence of the R2 ratio calculated with the asymptotic distribu- 
tion amplitude with (solid thick line) and without (dashed line) intrinsic 6-dependence. The 
thin line is as in Fig. 4. 
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